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We study the number statistics of ultracold bosons in optical Lattice using the slave particle 
technique and quantum Monte Carlo simulations. For homogeneous Bose-Hubbard model, we use 
the slave particle technique to obtain the number statistics near the superfluid to normal-liquid phase 
transition. The qualitatively behavior agree with the recent experiment probing number fluctuation 
[Phys. Rev. Lett. 96, 090401 (2006)]. We also perform quantum Monte Carlo simulations to ID 
system with external harmonic trap. The results qualitatively agree with the experiments. 
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I. INTRODUCTION 

The ultracold atoms in optical lattices have opened a 
new windows to investigate the strongly correlated sys- 
tems with highly tunable parameters^. The basic physics 
of these ultracold atoms is captured by the Bose Hubbard 
model, whose most fundamental feature is the existence 
of superfluid to Mott-insulator phase transition at zero 
temperatureSi^. In a very shallow optical lattice, the ul- 
tracold bosons are in superfluid phase and can be well de- 
scribed by a macroscopic wave function with long-range 
phase coherence^. In this case, the phase fluctuation is 
zero and the on-site number fluctuation is large. When 
the optical lattice is very deep, the bosons enter the Mott- 
insulating phase with flxed number of atoms per site and 
without phase coherence, i.e., the on-site number fluctu- 
ation is zero and the phase fluctuation is large"*'^. The 
physics of the MI phase is that, when the repulsive in- 
teraction between the atoms is large enough, the number 
fluctuation would become energetically unfavorable and 
the system would be in a number-squeezed state. This 
interaction induced MI phase plays an important role 
in the strongly correlated systems, as well as in various 
quantum information processing schemes?. 

In the past, some ultracold-atom experiments have 
been performed to detected these number-squeezed 
MI phase through the observation of increased phase 
fluctuations^'^'^ or through an increased time scale for 
phase diffusion^. Recently, the continuous suppression of 
on-site number fluctuations was directly observed by Fab- 
rice Gerbier et al by monitoring the suppression of spin- 
changing collisions across the superfluid/Mott-insulator 
transitionSi. By using a far off-resonant microwave field, 
the spin oscillations for doubly occupied sites can be 
tuned into resonance and the amplitude of spin oscilla- 
tion is directly related to the probability of finding atom 
pairs per lattice site. It was shown by Fabrice Gerbier et 
al that, for small atom number, the oscillation amplitude 
is increasingly suppressed with increasing lattice depths 
and completely vanishes for large lattice depths. In the 
MI region, this suppression persists up to some threshold 



atom number. The authors also compared their experi- 
mental results with the prediction of the Bose Hubbard 
model within a mean-fleld approximation at zero tem- 
perature. 

In this paper, we try to use the other approaches 
to study the number fluctuation beyond the zero- 
temperature mean-field theory. We first use the slave 
particle technique to obtained the number statistics at 
the critical points of the superfluid to normal liquid phase 
transition. The qualitative behaviors are the same as the 
recent experiment^ . In the second part of this paper, 
we step out the mean-field theory and perform quantum 
Monte Carlo simulation to the ID ultracold bosons with 
external harmonic trap. The numerical results reproduce 
the qualitative behaviors of the experiment. 

This paper was organized as follows. In Sec. ^1 we will 
describe the slave particle technique to the homogeneous 
Bose-Hubbard model. In Sec. IIIII we will perform quan- 
tum Monte Carlo simulation to the ID ultracold bosons 
with external harmonic trap. In Sec. IIVI wc will give our 
conclusions. 



II. SLAVE-PARTICLE APPROACH TO THE 
NUMBER FLUCTUATION OF HOMOGENOUS 
BOSE-HUBBARD MODEL 

We consider an ultracold atomic gas trapped in 
an three-dimensional optical lattice potential, Vo(r) = 
Vo^j^iSin^ {kvj), with wave vectors k — 2tt/X and A 
the laser wavelength. In real experiments, an additional 
harmonic potential is superimposed to the lattice poten- 
tial; however, we only pay attention to the homogeneous 
case in this section, which can be described by the fol- 
lowing homogeneous Bose-Hubbard Hamiltonian^, 

<ij> i i 

Here a| is the creation operator at site i, ni — aja; is 
the particle number operator, and {ij) denotes the sum 
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over nearest neighbor sites, t and U are the hopping 
amphtude and on-site interaction, respectively, 



t 

U 



J drw*{r ~ rO (^-^V^ + Voir)^ w{r - r^), 
9 J dr\w{r)\\ (2) 



In the following, we will use the slave particle technique 
to obtain the finite temperature number fluctuation at 
the critical points. In the slave particle languageiSiii, the 
bosonic creation operator a| and annihilation operator 
can be decomposed into 

aj = ^ \/a~fT\a + l)ii{a\, 

a=Q 

ttj = y^^Va + l\a)ii{a + 1|, (3) 

0=0 

where \a)i is the eigenstate of the particle number opera- 
tor Tii — alai with a the eigenvalue. In the slave particle 
language, every occupation state is identified as a type 
of slave particle, i.e., \a)i and i{a\ are mapped to a]^ ^ 
and aa,i, which are the slave particle creation and anni- 
hilation operators, respectively. Then al and can be 
rewritten as 



(4) 



Q=0 



The slave particle operators are defined to satisfy the an- 
ticommutation relation {aa,i,alj — SapSij in a slave 
fermion approach, and to satisfy the commutation rela- 
tion [aa,i,aj^j] = SafjSij in a slave boson approach. In 
order to reproduce the original bosonic commutation re- 



lation [ai,aj 



the slave particle operators must 



obey the constraint: 



(5) 



a=0 



a=0 



Substituting the slave particle transformation Q into 
Eq. Q, the Bose-Hubbard Hamiltonian can be replaced 
by 



<ij> a, (3 



(6) 



Following the steps in RefsiiSiH, we write the parti- 
tion function as an imaginary time coherent state path 
Integra! 
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(7) 



\ ( U 

S[aa,aa,X]^ / dr J 9^ - + — a(a - 1) 



t J2 E^/^+1^^ 



<ij> a,p 



(8) 



where Oq, i and introduced as ordinary complex 

numbers in the slave boson approach, and as Grassmann 
variables in the slave fermion approach satisfying the 
Grassmann algebra. The Lagrange multiplier field Ai(T) 
comes from the constraint ©, namely, Hi HJ2a ~ 
The unit has been set to ft = A;b = 1 in all formulas. 
In order to decouple the hopping term, we perform a 
Hubbard- Stratonovich transformation 



L'<i>*D<i>cxp 



- drt ($* - ^ + laa+i,iaa,t) 

<ij> a 



^ E + laa.jaa+i.j) 



(9) 



The Hubbard-Stratonovich field introduced here can 
be identified as the order parameter of superfluid for 
i^i) = (Y^a + laa,iaa+i,i)) = (oi). Wc then perform 
a Fourier transform on all the fields Ai by 



(10) 



where L is the total number of sites of the optical lattice 
and LOn is the Matsubara frequency, which equals (2n ~\- 
VjTr/ [3 or 2mT/ fi for fermionic or bosonic fields. After 
relaxing the constraint jS)) to one slave particle per site 
on average over the whole lattice, i.e., replacing Ak.n with 
a constant Ao.Oi we arrive at an effective action divided 
into two parts, 

5e//[$,aa,A] =5o-f 5/, (11) 
= iLfiX + X! E "k,n [-i^n + c(a)] ak,„ 

k.n a 

+ ^ek|<i>k,„|' = ^o'' + E^k|<i>k,„P, 



E E- 

k,k',Ti,n' a 



L/3 



(Va + la"k+k'), («+«')<") 



(Va + lak,„a"k+k'),(„+n')) 



*k' 



where c(a) = —iX—ai.j,+a{a—l)U/2, X = Xofi/ ^ L(3, and 
Ek = 2iX]f=i cos(fcia) with d and a being the dimension 
and spacing constant of the lattice. 
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Near the critical point, the order parameter $ is small 
and the perturbation can be performed in terms of Sj. 
The partition function Zq of non-interacting slave parti- 
cles comes from the contribution of the zeroth-order term 
Sq^ which is given by 



Zo 



-130.0 



Da^Daae 



(12) 



Here VLq is the zeroth-order thermodynamic potential, 



no = i\L ± ^ X! ^^(-^ ^ ^" 



-/3c(q) 



(13) 



where the +(— ) sign corresponds to the slave fermion 
(slave boson). After expanding e~^''ff up to second or- 
der of Si and integrating out the slave particle field, we 
obtain the new effective action to second order of the 
order parameter field, 

^i5,e//[$*,$] =/3r!o-^$;;,„G-i(k,iw„)$k,„. (14) 

k,n 



The Green's function G(k,iujn) is defined by 



-,a + l 



aU' 



(15) 



(16) 



where is the occupation number and equal to 

exp{f3[~i\ ~afi + a{a - l)U/2]} ± 1 ' 

in which the + and — sign correspond to slave fermion 
and slave boson, respectively. 

The saddle point approximation to the constraint field 
A means: dil/dX = 0; the particle number conservation 
condition requires —dn/d/j, = N. The mean-field ap- 
proximation means all the fluctuations coming from the 
Green's function would not be considered in the above 
two conditions. Then the following two mean-field equa- 
tions can be derived, 



Q=0 



= 1, 



Q=0 



an = — 



N 
L 



(17) 
(18) 



where the n = N/ L is the average particle density. 

Combining Eqs. H15|) - (|18|l . we can obtain the num- 
ber fluctuation at critical temperature of superfluid to 
normal liquid phase transition. We show the results in 
Figffl One can see that the qualitative behavior agree 
with the recent experiment probing number fluctuation 
[Phys. Rev. Lett. 96, 090401 (2006)]. 



III. NUMERICAL RESULTS: QUANTUM 
MONTE CARLO CALCULATION IN ONE 
DIMENSION 

Although mean field theory are applicable in higher 
dimensions, its application to ID is questionable due to 



- U/zt=4 






FIG. 1: The average probability of a site to be occupied 
by two bosons n^/n as a function of average particle density 
n, for different interactions U/zt. All the calculations are 
performed at the critical temperature of the phase transition. 



the large quantum fluctuations. In this section, we will 
come to the numerical calculation and focus our study 
on the Boson— Hubbard model (1) in a harmonic trapped 
potential in ID optical lattice. The trapped potential we 



(19) 



where L is the chain length and Vt = 0.02i. The method 
we use is quantum Monte Garlo (QMC) simulations using 
the stochastic series expansion techniquei^ii^. In the sim- 
ulation, we set the lattice is large enough to neglect the 
boundary effects and the inverse of temperature (3 — lOOt 
in order to reach the ground state properties. 

Fig[21show rial pioi — 1,2,3,4) v.s. the average boson 
density p, where Ua is the average boson density for a 
bosons of the system. From the figure we can see, n\ 
decreases with increasing average boson density p, while 
713, '^■4 increase with it. They all change monotonously 
while rt2 is nonmonotonous. When p is small, n\ is quite 
large, i.e., most of the sites are one particle occupied. In 
this region, there is a n = 1 Mott plateau in the middle 
of the chain(see figEJ- As p increases, n\ decreases while 
77.2 increases slowly. At some critical pci , ni suddenly in- 
creases fast. This is because when p > pd, a superfluid 
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find that ni develop concaves in the center of the trap. 
And these concaves become more and more deep which 




FIG. 2: (Color on line) ria/p as a function of boson density 
p for various U /t at Vt = O.Olt. The chain size is L = 100. 
na{a = 1, 2, 3, 4) are average bonson density for a bosons on 
each site. 
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FIG. 4: (Color on line) Wa/p as a function of boson density 
p for various U /t at Vt = O.Olf. The chain size is L = 100. 
na[a = 1, 2, 3, 4) are average bonson density for a bosons on 
each site. 



correspond with the decreasing of the total number of 
ni in Fig|21 n2 evolves nonmonotonously. It increase 
suddenly with the appearing of the superfluid region in 
the center of the trap and develop concave in the center of 
the trap with the appearing of the second Mott plateau. 
With the increasing of boson number, the rt3 and 714 are 
not important when the total number of bosons are not 
large. 



IV. CONCLUSION 



FIG. 3: (Color on line) Local particle density rn as a function 
of site i for various p a.t U = lOt. 



forms in the middle of the n — 1 Mott plateau in the 
middle of the chain. As p increase further, rt2 reaches 
its maximum at pc2 and then decreases with increasing 
p. This corresponds the formation of the n — 2 Mott 
plateau in the middle of the chain(also see Fig|2|). When 
U increase, both pd and pc2 increase. From the above 
discussion we can see that, 712 can be used to describe 
the Mott-superfluid transition. Our results are qualita- 
tively agree with Gerbier's experiment and the mean field 
theory above. 

In FigOand Fig^ we show the local particle density 
and onsite number statistics for various p at U — lOt. 
From Figl^l we can see that the formation of the n = 
1 and n — 2 plateau with increasing p. In Fig^l we 



We studied the number fluctuation of ultracold bosons 
in optical Lattice using the slave particle technique and 
quantum Monte Carlo simulation. By using the slave 
particle technique, we obtained the number statistic at 
the critical points of superfluid to normal liquid phase 
transition. We also performed QMC simulations to the 
one dimension Bose-Hubbard model with external har- 
monic trap. The numerical results qualitatively agree 
with the recent experiment^. 



Acknowledgments 

This work was supported in part by Chinese National 
Natural Science Foundation. The simulations were per- 
formed on the IIP-SC45 Sigma-X parallel computer of 
ITP and ICTS, CAS . 

Note adrferf.-Recently we became aware of a parallel 
numerical work^^ that reach similar conclusions. 



5 



^ I. Bloch, Physics World, 17, 25-29 (2004); I. Bloch, Nature 

Physics 1, 23-30 (2005). 
^ M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. 

Fisher, Phys. Rev. B 40, 546 (1989). 
^ D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner and P. 

Zoller, Phys. Rev. Lett. 81, 3108 (1998). 
" M. Greiner, O. Mandel, T. Esshnger, T. W. Hansch, and 

I. Bloch, Nature (London) 415, 39 (2002). 
^ C. Orzel, A. K. Tuchman, M. L. Fenselau, M. Yasuda, and 

M. A. Kasevich Science 291 2386 (2001). 
® P. Rabl, A. J. Daley, P. O. Fedichev, J. L Cirac, and P. 

Zoller, Phys. Rev. Lett. 91, 110403 (2003); G. Pupillo, A. 

M. Rey, G. K. Brennen, C. W. Clark, and C. J. Wilhams, 

J. Mod. Opt. 51, 2395 (2004). 

T. Stoferle, H. Moritz, C. Schorl, M. Kohl, and T. 
Esslinger, Phys. Rev. Lett. 92, 130403 (2004). 

* M. Greiner, O. Mandel, T. W. Hansch, and I. Bloch, Na- 
ture (London) 419, 51 (2002). 

® Fabrice Gerbier, Simon Foiling, Artur Widera, Olaf Man- 
del, and Immanuel Bloch, Phys. Rev. Lett. 96, 090401 



(2006). 

^'^ D. B. M. Dickerscheid, D. van Oosten, P. J. H. Denteneer, 
and H. T. C. Stoof, Phys. Rev. A 68, 043623 (2003). 

" Yue Yu and S. T. Chui, Phys. Rev. A 71, 033608 (2005); 
X. C. Lu, J. B. Li, and Y. Yu, Phys. Rev. A 73, 043607 
(2006). 

J. W. Negele and H. Orland, Quantum Many-Particle Sys- 
tems (Addison- Wesley, Redwood City, 1988). 
" H. T. C. Stoof, cond-mat/9910441 

A. W. Sandvik and J. Kurkijarvi, Phys. Rev. B 43, 5950 
(1991); A. W. Sandvik, J. Phys. A 25, 3667 (1992); A. W. 
Sandvik, R. R. P. Singh, and D. K. Campbell, Phys. Rev. 
B 56, 14510 (1997). 

^5 A. W. Sandvik, Phys. Rev. B 59 R14157701 (1999); O. F. 
Syljuasen and A. W. Sandvik, Phys. Rev. E 66, 046701 
(2002). 

B. Capogrosso-Sansone, E. Kozik, N. Prokofev, B. Svis- 
tunov, ,cond-mat /0609600, 



